Condensate density of interacting bosons: a functional renormalization group 

approach 
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We calculate the temperature dependent condensate density p°{T) of interacting bosons in three di- 
mensions using the functional renormalization group (FRG) . From the numerical solution of suitably 
truncated FRG flow equations for the irreducible vertices we obtain p°(T) for arbitrary tempera- 
tures. We carefully extrapolate our numerical results to the critical point and determine the order 
parameter exponent j3 ~ 0.32, in reasonable agreement with the expected value 0.345 associated 
with the XY-universality class. We also calculate the condensate density in two dimensions at zero 
temperature using a truncation of the FRG flow equations based on the derivative expansion in- 
cluding cubic and quartic terms in the expansion of the effective potential in powers of the density. 
As compared with the widely used quadratic approximation for the effective potential, the coupling 
constants associated with the cubic and quartic terms increase the result for the condensate density 
by a few percent. However, the cubic and quartic coupling constants flow to rather large values, 
which sheds some doubt on FRG calculations based on a low order polynomial approximation for 
the effective potential. 

PACS numbers: 05.10.Cc, 05.30.Jp, 03.75.Hh 



I. INTRODUCTION 

The condensed phase of interacting bosons is of- 
ten studied using Bogoliubov's celebrated mean-field 
approximation, 1 - However, if one tries to go beyond the 
Bogoliubov approximation and includes fluctuation cor- 
rections perturbatively, some terms in the perturbation 
series for the single-particle Green function diverge. The 
upper critical dimension for these divergencies is D = 3 
at zero temperature and D = 4 at finite temperature. 
In the critical dimension the divergencies are logarith- 
mic^ whereas in lower dimensions one encounters even 
stronger power law singularities. Physically, these diver- 
gencies arise due to the coupling of transverse fluctua- 
tions to longitudinal ones in the condensed phase; the 
gapless nature of the transverse fluctuations associated 
with the Goldstone modes gives then rise to singularities 
in the perturbation series of longitudinal correlation ianc- 
tionsj^SiSiZiS which have to be re-summed to all orders 
in perturbation theory to obtain meaningful results. The 
controlled calculation of physical properties of interacting 
bosons requires therefore non-perturbative methods. 

Recently several authors have studied inter- 
acting bosons by means of the renormalization 
group jVi 11 : 12 ! 13 ' 14 which is an efficient method to 
re-sum the perturbation series and remove the sin- 
gularities encountered in finite order perturbation 
theory. While most calculations so far have focused 
on properties of the superfluid ground-stat o 9 ' 10 i n 
or on the single-particle spectral function at zero 
temperature ) 12 i 13 Floerchinger and WetterichM have 
used the non-perturbative functional renormalization 
group (FRG) to calculate thermodynamic observables of 



the interacting Bosc gas in three dimensions. They used 
a truncation of the formally exact FRG flow equation for 
the generating functional T of the irreducible vertices^ 
based on the derivative expansion, ^ retaining field 
gradients and density fluctuations in the expansion of 
the effective potential up to second order. In Sec. [Til we 
shall present an alternative derivation of the resulting 
finite temperature flow equations based on the vertex ex- 
pansion*^ In contrast to Ref. [l4| , we shall then carefully 
analyze the temperature dependence of the condensate 
density p° (T) in the critical regime and extract the order 
parameter exponent /3 from the numerical solution of 
our truncated FRG flow equations. Our result (3 w 0.32 
is quite close to the expected value f3 « 0.345 of the 
three-dimensional XY-universality classJ^ We therefore 
conclude that our simple truncation of the vertex 
expansion yields quantitatively accurate results for the 
condensate density for all temperatures, including the 
critical regime. 

In order to estimate the effect of higher order many- 
body interactions, we shall in Sec. lIIII go back the deriva- 
tive expansion approac h 14 ' 16 and calculate the conden- 
sate density at vanishing temperature within a truncation 
to second order in the derivatives but to fourth order in 
the expansion of the effective potential U(p) in powers of 
density fluctuations p — p°, 



U(p) « 



4 rr(fc) 



fc=2 



(1) 



We find that the coupling constants and flow 
to rather large values, which sheds some doubt on the 
quantitative accuracy of calculations based on the gradi- 
ent expansion with quadratic approximation for the ef- 
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fective potential. Nevertheless, we find that the three- 
and four-body interactions described by and 
have only a rather small effect on the numerical value of 
the condensate density. 



so that the infrared cutoff A suppresses long-wavelength 
fluctuations and we recover our original model for A — > 0. 
For convenience we use the Litim regulator— 

R A (k) = (1 - 5 kfl ) Zl 1 (e A - e k ) 6 (A 2 - k 2 ) , (6) 



II. FRG APPROACH TO THE CONDENSED 
BOSE GAS 

A. Vertex expansion of the FRG flow equation 

Starting point of our investigation is the following Eu- 
clidean action describing bosons with mass m subject to 
a repulsive contact interaction uq, 



S$,il>] 



dr 



i)(r, T)(d T - - v)4>(r, t) 



(2) 



where the chemical potential /i and the inverse temper- 
ature (3 ~ 1/T are fixed. The spatial integrals should 
be regularized by means of a short-distance cutoff A^" 1 , 
which is related to the finite extent of the interaction or, 
for hard core bosons, to the size of the particles. The 
model (J2J) depends on three dimensionless parameters 
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where at this point we do not specify the dimensionality 
D of the system, and the factor of 2tt in the definition of 
the dimensionless temperature T is introduced for later 
convenience. We focus on the condensed phase, where 
the global [/(l)-symmetry of the action ^ is sponta- 
neously broken and the field ip has a finite expectation 
value (jP = (ip(r, r)), which by translational invariance is 
independent of space r and imaginary time r. Without 
loss of generality, we choose <fi° to be real. 

To derive formally exact FRG flow equations for the 
one-particle irreducible vertices of our model, we add a 
cutoff dependent regulator function i? A (fc) to the inverse 
free propagator in the Gaussian part of the action © . In 
momentum-frequency space the inverse free propagator 
is then 



(4) 



where ito is a bosonic Matsubara frequency, e k — k 2 /2m 
is the free dispersion in momentum space, and K — 
(k, iui) is a collective label. The regulator function i?A(&) 
should satisfy 
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where the dimensionless wave-function renormalization 
factor Za is defined in Eq. ([23]) below. The cutoff depen- 
dent irreducible vertices T^' m ' (K[, . . . , K' n ; K m , . . . , Ki) 
in the condensed phase are defined via the func- 
tional Taylor expansion of the corresponding generat- 
ing functional^ in powers of the fluctuations 
<5>k - 5k.o4>°, 
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The derivative of the functional T\[4>, cjj] with respect to 
the infrared cutoff A can be expressed in closed form in 
terms of a deceptively simple FRG flow equation ) 15 i 16 i 17 
which is equivalent to an infinite hierarchy of integro- 
differential equations for the cutoff-dependent vertices 



^(n,m) 



K' ■ K 



,ifi). Following Refs. [21.2 



23[, we fix the flowing order parameter (j) A = {ip(r,r)) 
by demanding that the vertices r^ 1 ' - 1 and T^' 1 ) w ith a 
single external leg should vanish identically for any value 
of the cutoff A. 

We are interested in the temperature dependent order 
parameter <f>1 = limA^o <Pa , which determines the con- 
densate density via p° = (0°) 2 . The exact FRG flow 
equation for the flowing order parameter <j) A depends on 
the flowing normal and anomalous self-energiesj2I 



T k 1a \k,K) = I%(K) 



Tf 2 \K 1 -K) = Tf'»{-K 1 K) = I%(K), 



_ r (2,o) 



(8) 
(9) 



and on the four types of vertices with three external legs, 

r (3,0) 5 pp,l) ) p(l,2) and p(0,3)_ Tq calculate the order pa _ 

rameter we therefore need the flowing self-energies and 
the flowing three-legged vertices, whose flow equations 
depend again on higher order vertices with four and more 
external legs. To obtain a closed system of FRG flow 
equations, we shall use here the truncation proposed in 
Ref. [l2j], which amounts to the following parameteriza- 
tion of the non-zero vertices with three and four external 
legs, 



T^' 2> (K[, K' 2 ; K 2 , Ki) 



4>l[u A {K 1 ) + u h {K 2 )], (11) 

UAiKi-KJ+UAiK'z-K!). 

(12) 
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The system of flow equations for the order-parameter 
and the self-energies is closed by demanding that the 
momentum- and frequency dependent function u A (K) is 
related to the flowing self-energies via 



(K) 



a A (K) + p° A [u A (0) 
P>a(K), 



u A (K)], (13a) 
(13b) 



where "' 1 



Pa = (<^a) 2 i s the flowing condensate density 
and a A (K) is another if -dependent function satisfying 
<7 A (0) =li- p A u A (0). Eqs. (fTOl - [T2"j) relate the vertices 
with three and four external legs to the normal and 
anomalous components of the irreducible self-energy, and 
thus close the FRG flow equations for the order parame- 
ter and the self-energies. Our parameterization of the 

IT3bll 



self-energies given in Eqs. (|13al I13b[) is motivated by 
the pioneering insights gained by Nepomnyashchy and 
Nepomnyashchy, 2 implying that only the contribution to 
the self-energy contained in the function u A (K) exhibits 
a non-analytic if-dependence for A — > 0, while the con- 
tribution o A (K) remains analytic. As noted in Ref. 
the above truncation satisfies the Hugenholtz-Pines rela- 



tion 
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££(0)-E£(0)= M , 



as well as the Nepomnyashchy identity, 



lim Ef (0) = 0, 



(14) 



(15) 



which holds at finite temperature for D < 4 and at zero 
temperature for D < 3, since in these cases 



lim u A (0) = 0. 



(16) 



The truncation (fTOHT^I amounts to the following approx- 
imation for the generating functional T A [<p, 4>] defined in 
Eq. 0, 



p K u A (K)p. 



K ■ 



r A [& 0] - Tf + I j> K a A (K)cjy K -r - 

.Ik .1 /, 

(17) 

where px — Jq 4>q4>q+k are the Fourier components of 

the density^ and T A is an interaction correction to the 
grand canonical potential in units of the temperature. 
As shown in Ref. 12], for A — > the function u A (K) 
develops a non-analytic dependence on K, so that the 
corresponding effective potential in real space and imag- 
inary time is non-local. 

To further simplify the FRG flow equations, we shall 
replace on the right-hand sides of the flow equations, 



u A (K) -> u A (0) = u A . 



(18) 



Within this truncation, the FRG flow equation for the 
condensate density p A — ((f> A ) 2 reduces to 



OaP° a 



[2G N A (Q) + Gf{Q)], 



(19) 



while the normal and anomalous components of the self- 
energy satisfy 12 



d A ^(K) = 2u A j [G N A (Q) + G£(Q)} - 4u 2 A p° A [G N A (Q) [G%(Q + K) + G N A {Q-K) + G N A (-Q + K) 

+2G£(Q - K)] + 2Gi(Q) [G A (Q + K) + G a i (Q + K)] } , (20) 
d A Y,i{K) = 2u A f G N A (Q)-4u\pl f {G^{Q)[G N A (Q + K) + G N A (Q-K) + Gi(Q + K) + Gi(Q-K)\ 

JQ J Q L 

+6£(Q) [G% (Q + K) + G% (Q-K) + 3G£(Q + K)]}. (21) 



Here the single-scale propagators G A (K) and G A {K) are where G A 1 (if) = G \ {K) — H A (K), and 
defined via the matrix equation 



G»{K) G£(K) 

g£(k)* g%(-k) 



= -G A {K)[d A G Q \{K)]G A {K), 

(22) 



G 0>A (K) = 
V A (K) = 



G 0<A (K) 
G 0iA (-K) 



(23) 
(24) 
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B. Low-energy truncation and results 

Following Ref. [l^] we expand the analytic part a A (K) 
of the self-energy in powers of momenta and frequencies 
up to quadratic order, 



a A (K) a n(l-X A )+iu(l-Y A ) 
+ e k {Zl x -l)-(iw)V A , 



(25) 

where X A = p A u A / p,. Note that at the initial scale 
X Ao = *Aq = z a = 1 and Vaq = 0, so that a Ao (K) 
vanishes. While at zero temperature it is essential to 
retain the couplings Y A and V A associated with the fre- 
quency dependence of the self-energy, these couplings are 
irrelevant at the critical fixed point associated with Bosc- 
Einstein condensation, which is a classical phase transi- 
tion. Since in this section we are interested in the finite 
temperature behavior of the condensate density, for our 
purpose it is sufficient to retain only the couplings X A 
and Z A in Eq. (j2"5j) . 

a A {K) « <r A (fe, iuj = 0) « (i - r A + e h {Z A l - 1), (26) 

where we have introduced the notation 



rA 



p° A u A . 



(27) 



With this truncation, the normal and anomalous propa- 
gators are simply 



Gf{K) 
Gi(K) 



Z A 1 € k 



r A + R A (k) 



[Z. 



A e * " 

r A 



r A -i? A (fe)] 2 



^ -rl + [Z]r l e k + r A - R^k)] 2 ' 
while the corresponding single-scale propagators are 
G%(K) = ~[d A R A (k)] 



(28) 
(29) 



[~iuj - Z A 



(k-r A 



RA{k)f 



[Zl\ k + r A -R A {k)Y\ 



G A {K) = -[d A R A {k)\ 



2r A [R A (k)-Z] v 1 e k -r A } 



(30) 



(31) 



[u 2 -r\ + [Z A 1 e k + r A -R A {k)] 2 \ 

Substituting these expressions into Eqs. (fT9l[2"Tj) and us- 
ing the Litim regulator (O, the momentum integrations 
and frequency summations appearing in the FRG flow 
equations for the three independent couplings p\ , u A and 
Z A at finite temperature can all be performed analyti- 
cally. In D dimensions the final result can be written in 
the following form, 

D+2 



A9 A/ 5 A = 



2Kj 



A_ 

Ao 



irD 

3 S o , 2 (PE A )P 1 



1 - 



VA 



Ad A u A 



8K D 
' ttD 



(2) 
A 

D+2 



D + 2 

PS 2 , 2 0E A 

VA 



(32) 



D + 2 



~2 rr- 

u A Z A 



x pS , 3 ([3E A )P { A 3) - ^S 2 , 3 ((3E A )P A 



(i) 



(33) 



Ad A Z A = r] A Z A , 
where the flowing anomalous dimension is 



(34) 



VA 



D+2 



^(A) / v; ( -/ A z; l .r.s,oi ie \ i. { :>,+ 



Here K D = 2 1 - d tt- d / 2 /T[D /2] is the surface of the D- 
dimensional unit sphere divided by (2ir) D , the dimension- 
less inverse temperature (3 = l/T is defined via Eq. (|3b[) . 
and we have introduced dimcnsionlcss quantities 



Pa = P>o°, 

u A = 2itlu a Aq~ 2 , 

f A = p A u A = 2mp° A u A /Al, 

E A = \/e A (i A + f A ), 

~e A = Zl l A 2 /Al 

as well as dimensionless coefficients, 



P A ] = 7^ + ^eA, 



7. 
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,(2) 



f\ + r A e A + e A , 



(36a) 
(36b) 
(36c) 

(36d) 
(36e) 



(37a) 
(37b) 



P 



(3) 



r A + ~r A £A + \f A ei + -e\. (37c) 



Finally, the dimensionless functions S kt i(x) are defined 
in terms of the bosonic Matsubara sums, 



Sk,i{x) 



oo 

E 



(n 2 



„2\l ■ 



(38) 



which can be expressed in terms of the Bose function and 
its derivatives. 

The system of first order differential equations given 
by Eqs. (|32n34"|) can easily be solved numerically. It 
is convenient to consider quantities as functions of the 
logarithmic flow parameter I = — ln(A/Ao), renaming 
P° Age -i — * p\, and analogously for the other couplings. In 
Fig. [1] we show the typical RG flow of the couplings {i/ 
and the flowing anomalous dimension r\i for three differ- 
ent temperatures (above, below and at the critical tem- 
perature). As expected, only at the critical temperature 
the flowing anomalous dimension has a finite limit (which 
can be identified with the critical exponent 77), while for 
T < T c the condensate density approaches a non-zero 
value p®. The temperature dependence of the corre- 
sponding order parameter (f>® = y^52 is shown in Fig. [2] 
for three different values of the interaction strength. In 
order to extract the order parameter exponent (3 from 
our numerical results shown in Fig. [2l one should care- 
fully fit the curves in a sufficiently small temperature 
interval below the critical temperature to a power law, 



a (f c - fy 



(39) 



Due to the lack of a priori knowledge about the proper 
interval for the power-law fit, and due to the finite ac- 
curacy of the numerical data, it is non-trivial to extract 
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f = 


43 


f = 


45.1318 


f = 


47 




FIG. 1: (Color online) Typical FRG flow of the condensate 
density, the anomalous dimension and the dimensionless in- 
teraction ui for D = 3 and p, = Ho — 1. Solid line: criti- 
cal temperature (T c = 45.1318 for our choice of parameters); 
dashed line: f = 47 > f c \ dotted line: f = 43 < f c . 



the critical exponent (3 from the numerical solutions of 
the FRG flow equations shown in Fig. [TJ In fact, in a 
recent FRG calculation^ of the temperature dependent 
condensate density the critical exponent j3 was not deter- 
mined, apparently due to a lack of numerical accuracy. 
Here we present an extrapolation procedure which allows 
us to obtain the critical exponent (3 with high accuracy. 
The crucial point is that one should use a series of in- 
creasingly narrow intervals close to the critical point for 
the fitting procedure. Specifically, we use intervals of the 
form 



I z = [(l-2- 2 )f c ,f c ], 



(40) 



which are parameterized in terms of the zoom factor z. 
For increasing values of z we fit the data in the corre- 
sponding temperature interval I z to the power law (139|) 
and extract (3(z). This procedure is illustrated in Fig. [3] 




FIG. 2: (Color online) Temperature dependence of the or- 
der parameter (j>i = y/32 m three dimensions for jl = 1 and 
three different values of the dimensionless bare interaction ito 
defined in Eq. (|3c)) . 



We then plot our results for j3(z) as a function of z and 
extrapolate for z — > oo. As shown in Fig. [3l the depen- 
dence of P(z) on the zoom factor z can be described by 



(3(z) w (3[1 - e 



-a(z+7)i 



(41) 



with some non-universal numbers a and 7. The critical 
exponent /3 can then be identified with {3 = lim 2 ^oo j3(z). 
Given the simplicity of our truncation our final result 
(3 fss 0.32 is in reasonable agreement with the accepted 
value (3 « 0.345 of the three-dimensional XY-universality 
class. 18 We have checked that our extrapolated result for 
(3 is independent of the various non-universal parameters 
of our model, such as Ao or the value of uq. We therefore 
believe that, in spite of its simplicity, our truncation of 
the vertex expansion yields accurate results for the con- 
densate density for all temperatures. 



III. DERIVATIVE EXPANSION WITH 
QUARTIC EFFECTIVE POTENTIAL 

The truncation of the vertex expansion used in the 
previous section is equivalent to approximating the gen- 
erating functional T\[(f>, <p] of the irreducible vertices by 
Eq. ([17]) . A weak point of this truncation is that it ne- 
glects many-body interactions involving more than two 
particles, encoded in the irreducible vertices with more 
than four external legs. Within the framework of the ver- 
tex expansion it is rather difficult to take these vertices 
into account. In this section we shall therefore use the 
derivative expansion 16 with cubic and quartic terms in 
the effective potential to estimate the effect of higher or- 
der many-body interactions. For simplicity we consider 
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1 - 2- 2 )T ( 



0.32 r 



0.31 



0.3 




fit function: 
P = 0.3209 



( 




(1 - 2-^ +1 ))f c f c T 

FIG. 3: (Color online) Iterative procedure to extract the or- 
der parameter exponent from the FRG results. We fit our 
numerical FRG results in a series of intervals I z given by 
Eq. (|40p to a power law (I39|l and thus determine the critical 
exponent P(z) for a given value of the zoom factor z. The z- 
dependence of P{z) is then extrapolated for z — * oo as shown 
in Fig. H 



0.29 1 1 1 1 

5 10 15 20 

zoom factor z 

FIG. 4: (Color online) The dependence of P(z) on the zoom 
factor z can be described by the indicated fit function [see 
also Eq. (@TJ] with a = 0.1945 and 7 = 8.545. The data are 
based on the numerical solution of the FRG flow equations 
for iio = p, = 1. 



to write the local effective potential in the form 

Ua( P ) = C/f + f (P-Plf. (44) 

Here 



U^^-^iplf, (45) 



in this section only the quantum renormalization of the 
condensate density at vanishing temperature. 

If we approximate u A (K) » ua(0) = ma in Eq. (fTT)) 
and use the low-energy expansion (|25|) for the analytic 
part <j{K) of the self-energy, then our truncated vertex 
expansion of Sec. Ill Bl amounts to the following approx- 
imation for the generating functional of the irreducible 
vertices, 



T A [<l>,<l>]-ii I d D r I drp(r,r) 

U A (p(r, r)) + (1 - Y A )ftr, r)5 T 0(r, r) 



d D r I dr 



<z A 



\V6(r, 



2m 



V A \d T( i>(r,T)\ 



(42) 



where p(r,r) = \<j>{r , t)| 2 , and we have used the identity 
Kl~X A )p+ ^ = w+fip-plf - f(plf, (43) 



where V is the volume of the system. Recall that through- 
out this work we normalize Fa[<^,</>] such that it re- 
duces to the two-body part of the bare action at the 
initial RG scaled where X Ao = Y Ao = Z Ao = 1 and 



V Ao 



_ r(o) _ 



(the chemical potential term is included 



in the non-interacting Green function). 

To investigate the effect of higher order many-body in- 
teractions, we now generalize the above ansatz by replac- 
ing the effective potential by a fourth order polynomial 
in the density 



U A (p) 



u 



(2) 



2! 



U 



(3) 



U 



(4) 



3! 



{p-pir + ^ip-pir- (46) 



In fact, for the truncated generating functional Fa[<^, <f>] 
of the form (|4"2"1) with an arbitrary local effective po- 
tential U A (p) the exact FRG flow equatio n 15 ' 16 implies 
the following partial differential equation for the effec- 
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tive potential 



14 



djj 



A D+1 



D + 2J 2mZ A 
P U"{p) + U'(p) + V A u 2 



2,,, 2 ' 



2tt [pU"(p) + U'{p) + V a lu 2 } 2 - [pU"{p)] 2 + Y 2 uj 



(47) 

At finite temperature, the frequency integral should be 
replaced by a bosonic Matsubara sum, J |^ — » T 1 ^^. To 
obtain an approximate solution of the partial differential 
equation (14"?| . we expand U A (p) in powers of p — p A . The 
flowing condensate density p A is then determined by 



5C/ A (p) 



dp 



and the expansion coefficients are 



u A 0) = uM) 



u 



(fc) _ d fc C/ A (p) 



(48) 



(49) 



Pa 



Taking derivatives of Eq. |47|) , we obtain the flow equa- 
tions for the condensate density and expansion coeffi- 
cients, 



dAP° A 



-^§- p (dKU A {p))\ pV 



'A 



dp 



-(d A u A ( P ))- 



u 



(fc+i) 



8 



C/i 2) d P 



{d A U A { P )) 



(50) 



(51) 



The flow of the couplings Z A ,Y A and Va. related to the 
single-particle Green function can be derived by insert- 
ing our ansatz ([42]) into the exact FRG flow equation 
for r A [0, (f>\ and comparing terms with the same num- 
ber of gradients on both sidesj 15 i 16 If we retain only the 

(2) 

quadratic coupling U A =u A and set Y A = 1 and V A = 
we recover the flow equations (|32H55)) obtained in Sec. [TT] 
within the vertex expansion. If we retain in addition 
the cubic and quartic couplings and U A ^ the result- 
ing system of equations is lather lengthy. In this work 
we do not give these equations explicitly because their 
derivation is straightforward and they can only be anal- 
ysed numerically anyway; technical details can be found 
in Ref. [26|. Note that the couplings U A and V A are not 
dimensionless; for our numerical analysis it is convenient 
to work with the corresponding dimensionless couplings 



U 



(fe) _ 



2mUT>K t 



(k) A D(k-l)-2 



V A = 



A 2 Va 
2m 



(52) 



For simplicity, we have explicitly solved the coupled 
flow equations for the seven dimensionless couplings 
U A \U A \U A \p A ,Z Al Y A and V A only at zero temper- 
ature. The typical RG flow for the coupling constants 
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FIG. 5: (Color online) RG flow of the dimensionless couplings 

U A 2 \U < £ ) and defined in Eq. J52J for m = 3, Jx = 1, 

D = 2 and T = 0. 



U/f* , U A ^ and in two dimensions is shown in Fig. [5] 
Obviously, during the RG flow the cubic and quartic cou- 
plings ujp and JJ 1 ^ become orders of magnitude larger 

~ (2) 

than the quadratic coupling U A . Keeping in mind that 
a low order truncated polynomial approximation for the 
effective potential is only justified if the neglected higher 
order terms are smaller than the retained lower order 
terms, we conclude from Fig. [5] that results based on 
truncations of the derivative expansion with quadratic 
effective potential should be considered with some skep- 
ticism. This is also the case in three dimensions where a 
similar calculation (not shown here) leads to a divergence 
of both coefficients U A and U A for A — > 0. 

In spite of this, the widely used quadratic approxima- 
tion for the effective potential can still give acceptable 
results for some physical quantities. For example, the 



r(4) 



large values of U A and U A shown in Fig. [5] lead only to 
a small correction for the condensate density, as shown in 
Fig. O Interestingly, the inclusion of cubic and quartic 
terms in the expansion of the effective potential renor- 
malize the condensate density to larger values, so that 
the quadratic approximation overestimates the effect of 
fluctuations. It would be interesting to determine the 
true form of the effective potential in two and three di- 
mensions by directly solving the partial differential equa- 
tion (T4T|) . This seems to be a rather difficult numerical 
task which is beyond the scope of this work. 



IV. SUMMARY AND OUTLOOK 

In summary, we have used two different truncation 
strategies for the formally exact FRG flow equation for 
the generating functional of the irreducible vertices to 
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FIG. 6: (Color online) RG flow of the condensate density for 
Uo = 3, p, = 1, D = 2 and T = with and without the cubic 
and quartic couplings and 

calculate the condensate density of the interacting Bose 
gas. Our first strategy presented in Sec.HJJ.is based on the 
truncated vertex expansion recently proposed in Ref . [l2[ • 
We have further simplified this truncation at finite tem- 
perature to obtain a closed system of flow equations for 
the condensate density, the effective interaction, and the 
wave-function renormalization factor at finite tempera- 
ture. These flow equations are equivalent to the flow 
equations recently derived by Floerchinger and Wetterich 
within the derivative expansion^ From the numerical so- 
lution of these flow equations we have obtained a quan- 
titatively accurate description of the critical regime. We 
have also developed an extrapolation procedure to ex- 
tract the order parameter exponent (3 with high accuracy 
and moderate computational effort from the numerical 
solution of the FRG flow equations. 

In order to estimate the renormalization of the conden- 
sate density by three-body and four-body interactions 



which are not included in our truncated vertex expan- 
sion of Sec. UH we have used in Sec. [Till the truncated 
derivative expansion in combination with a polynomial 
approximation for the effective potential and derived ap- 
propriate FRG flow equations at zero temperature. We 
have shown that the RG flow drives the cubic and quartic 
coefficients llf^ and in the expansion of the effec- 
tive potential U\(p) to rather large values. This indicates 
that a low order polynomial approximation is not appro- 
priate to describe the true form of the effective potential 
in the condensed phase. Although the inclusion of the 
coefficients and leads only to a small positive 
correction of the condensate density, other physical ob- 
servables might be more strongly affected by the three- 
and four-body interactions described by these couplings. 

By using both the vertex expansion and the deriva- 
tive expansion, we have clearly established the relation 
between these approximation strategies and have illus- 
trated their advantages. In this work we have focused on 
the condensate density, because it appears naturally as 
one of the flowing couplings in the FRG flow equations. 
Our work can be extended in several directions: First of 
all, by keeping track of the field-independent part 
of the generating functional in Eq. (|f 7p (which can be 
identified with the interaction correction to the grand 
canonical potential per unit volume) one can also ob- 
tain the FRG flow of any thermodynamic observable of 
interest^ In the critical regime our numerical extrapola- 
tion procedure outlined in Sec. UTI should again be useful 
to obtain quantitatively accurate results for other criti- 
cal exponents. It would also be interesting to generalize 
the vertex expansion approach developed in Ref. [f2j to 
finite temperatures and calculate the single-particle spec- 
tral function in the critical regime. 
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